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Abstract 

In this paper we explore stochastical and statistical properties of so-called recurring spike 
induced Kasner sequences. Such sequences arise in recurring spike formation, which is needed 
together with the more familiar BKL scenario to yield a complete description of generic space- 
like singularities. In particular we derive a probability distribution for recurring spike induced 
Kasner sequences, complementing similar available BKL results, which makes comparisons 
possible. As examples of applications, we derive results for so-called large and small curvature 
phases and the Hubble-normalized Weyl scalar. 
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1 Introduction 

It is no understatement to say that the work of Belinskh, Khalatnikov and Lifshitz (BKL) [TJ ® [3] 
has set the stage for much of subsequent investigations about the detailed nature of generic spacelike 
singularities. The starting point for their analysis is their 'locality' conjecture, which states that 
asymptotically toward a generic spacelike singularity in inhomogeneous cosmology the dynamics is 
local, in the sense that each spatial point is assumed to evolve toward the singularity individually 
and independently of its neighbors as a spatially homogeneous model [1] [3] . A second important 
conjecture of BKL is that some sources, like perfect fluids with sufficiently soft equations of state 
such as dust or radiation, lead to models that generically are asymptotically 1 vacuum dominated', 
i.e., asymptotically toward a generic spacelike singularity the spacetime geometry is not influenced 
by the matter content, even though, e.g., the energy density blow up [TJ[51[3]. This leads to a 
picture of temporally oscillating Kasner states (i.e., vacuum Bianchi type I solutions, which have 
flat spatial curvature) where two subsequent Kasner states are determined by the vacuum Bianchi 
type II solutions, so-called single curvature transitions in the nomenclature of [4]0 

Several recent papers @] — [ID] have provided further evidence for the BKL scenario, but it 
has also been found that there are timelines whose evolution is different, exhibiting 'non-local' 
recurring spike formation. Remarkably, BKL and recurring spike behavior are intimately linked in 
a hierarchical manner: by hierarchies of subsets and by a solution generating algorithm that links 
all building blocks for understanding asymptotics at generic spacelike singularities to the bottom 
of the hierarchy, which consists of the Kasner solutions. 

In a recent paper [5] it was shown that recurring spike formation could be analytically described 
by means of combining the exact inhomogeneous solutions found by Lim [7] into concatenation 
blocks that connect families of timelines affected by spike formation at common initial and final 
Kasner states. It was also shown that the concatenation blocks come in two versions, so-called 
high velocity transitions, TU\, and joint low/high velocity spike transitions, 7j , where the word 
transition, as in the BKL case, signifies a change from one Kasner state to another. Hence both 
BKL dynamics and spike dynamics can be described in terms of temporally oscillating Kasner 
states, and these oscillations can be described by discrete maps. 

The outline of the paper is as follows. In the next section we describe the BKL and spike maps, 
and the maps they in turn give rise to, and some of these maps properties; due to the central feature 
that all these maps describe changes in Kasner states, we refer to them collectively as Kasner maps. 
Furthermore, we derive a distribution function that allows one to statistically analyze properties 
associated with spike induced Kasner sequences, which complements previous corresponding BKL 
results. In Section [3] we explore some statistical features the maps give rise to. In particular we 
prove that statistically so-called small curvature phases dominate over large curvature phases, and 
we also give stochastical relations for the 'Hubble-normalized' Weyl scalar, for generic BKL as 
well as spike induced Kasner sequences. We conclude the paper with a few remarks about future 
possible research projects concerning explorations of chaotic features of recurring spike formation, 
which we leave for the interested reader. 

2 Kasner maps 

A Kasner state can be described in terms of the gauge invariant Kasner parameter u, see e.g. [2l|4l 
[TOl [TT] If] Next we describe the two types of discrete maps which determine the change in u due to 
the transitions, beginning with the map that arises from the Bianchi type II solutions (i.e., single 
curvature transitions) in the BKL scenario. Throughout, the time direction is chosen to be toward 
the initial singularity. 

The term transition denotes a change between two different Kasner states; 'single curvature' refers to a single 
degree of freedom that describes the spatial curvature, which is excited during the transition. 

2 The 1-parameter family of Kasner solutions is often given in terms of the line element ds 2 = — dt 2 + t 2pi dx 2 + 
t 2 ^dy 2 +t 2 P3dz 2 , where Pi + Pi + P3 = Pi + p\ + p\ = 1- The Kasner parameter u can be defined via P1P2P3 = 
-u 2 (l + u) 2 /(l + u + u 2 ) 3 . 
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2.1 The BKL map 

The single curvature transitions induce the following BKL map [21 lllj : 

if u 1 6 [2, oo ) , 




(1) 

if u 1 G [1,2] . 



This map is often referred to as the Kasner map, but since we here will use the term Kasner 
map as an umbrella term for maps that induce changes in a Kasner state, we instead refer to this 
particular Kasner map as the BKL map. 

2.2 The BKL era map 

Iterations of the BKL map generate, from every initial value uq 6 [1, oo), a (finite or infinite) BKL 
sequence of Kasner epochs (u;);=o,i,2,...- We say that the Kasner epochs ui and ui+\ belong to 
the same (BKL) era if = ui — 1, thereby leading to a partition of a BKL sequence of Kasner 
epochs into eras [2] [11], e.g., 

1.14 7.29 6.29 5.29 4.29 3.29 2.29 1.29 3.45 2.45 1.45 2.24 1.24 -> . 



Let us denote the initial (= maximal) value of the Kasner parameter u in era number s by u s , 
where s = 0, 1, 2, . . . . Following [2] [IT] , we decompose u s into its integer part k s = [u s ] and its 
fractional part {u s }, i.e., 

u s = fc a + {uj . (2) 

The number k s represents the (discrete) length of era s, which is simply the number of Kasner 
epochs it contains. The final (= minimal) value of the Kasner parameter in era s is given by 
1 + {u s }, which implies that era number (s + 1) begins with 

U * +1 = {£} • (3) 

The map u s i— >• u s+ i is (a variant of) the BKL era map, which we usually, for brevity, refer to as 
just the era map; starting from u = u it recursively determines u s , s = 0, 1, 2, . . ., and thereby 
the complete BKL Kasner sequence (ttj)/=o,i,... 01 epochs. 

The era map admits a straightforward interpretation in terms of continued fractions. Consider 
the continued fraction representation of the initial value, i.e., 

uo = k H = [k ;ki,k2,k 3 ,...] . (4) 

k 2 H 

The fractional part of uo is [0; k\, k 2l k^, . . .]; since Ui is the reciprocal of {uo}, it follows that 

ui = [ki;k 2 ,k 3 ,k 4 ,...} . (5) 

Therefore, the era map is simply a shift to the left in the continued fraction expansion, 

u s = [k s ; k s+ i, k s+2 , ■ ■ ■ ] t-t u s+i = [k s +i] k s+ 2, k s+3 , . . . ] . (6) 

In terms of continued fractions, the Kasner sequence (it;)z<=R generated by uo = [ko] fei, k-2, ■ ■ •] 
is given by 

u a = u 





h,k 2 , ■■■ 


] -► [ko 


- l;k 1 ,k 2l ... 


] ^ [ko 


-2;k 1 ,k 2 ,-. 


.]-.... 


[i;h,k2 


[fci; 


k 2 ,k 3 , ... 


\^[h 


- l;k 2 ,k 3 ,... 


]"►[*! 


~2;k 2 ,k 3 ,.. 


.]^... 


-> [i;^2,^3 


[k 2 ; 


k 3 , k^,... 


1 ^ 


- l; ^3, ki, . . . 


H[fe 2 


- 2; fc 3 , A: 4 , . . 


.]-..-. 





u 2 

Since u s G [fc s , fc s + 1), the number A: s represents the (discrete) length of era s, which is simply 
the number of Kasner epochs it contains. Therefore, passing on to the stochastical interpretation of 
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(generic) Kasner sequences of epochs, we find that the probability that a randomly chosen era s of a 
Kasner sequence (it/)i=o,i,2,... of epochs has length to £ N corresponds to the probability that k s = 
m, or, equivalently, to the probability that u s S [m, m+1). Since the sequence (ho, ki, fo, • ■ ■ ) arises 
as the continued fraction expansion of uq this probability in turn corresponds to the probability 
that a randomly chosen partial quotient in the continued fraction expansion is equal to m. Hence 
we resort to Khinchin's law [12], which states that the partial quotients of the continued fraction 
representation of a generic real number are distributed like a random variable whose probability 
distribution is given by 

*<»>-*> (^)-'-(^)- m 

which leads to 

Probability (length of era = n) = L{n) — K{n) . (8) 

It follows that 42% of the BKL eras of a generic BKL Kasner sequence of epochs contain merely 
one Kasner epoch; 17% of the eras contain 2 Kasner epochs, and, e.g., 1.4 * 10 _2 % of the eras 
consist of 100 Kasner epochs. 



2.3 The spike map 

Let us now consider the Kasner map that is induced by the spike transitions Th\ and 7] - It is 
shown in [5] that, remarkably, the two types of transitions give rise to the same Kasner map which 
we refer to as the spike map. This map is obtained by applying the BKL map twice, see [5] , which 
results in the following relations: 

u- — 2 w_ <E [3, oo) , 

. (U--2)" 1 u-€[2,3], 
" H " M^-l)- 1 -!)" 1 u-e [3/2,2], W 
(u--!)- 1 -! U- 6 [1,3/2]. 



2.4 The spike era map 

Asymptotically we expect that spike dynamics is reduced to dynamics along timelines associated 
with a spatial 'spike' surface. Furthermore, this gives rise to recurring spikes described by sequences 
of Th\ and 7j , which yield iterations of the spike map 9 . Iterations of the spike map ([9]) generate, 
from every initial value Uq € [l,oo), a (finite or infinite) recurring spike-generated sequence of 
Kasner epochs (u;);=o,i,2,...- We make a partition of (wi)z=o,i.2,... into recurring spike-induced 
eras, which we denote as spike eras, or, for brevity, eras: ui and it;+i belong to the same era if 
ui + i = ui — 2. If ui + i arises from ui by one of the other three laws of ©, we speak of a change of 
era. 




era era era era era 



We denote the initial (= maximal) value of the Kasner parameter u in era number s (where 
s = 0, 1, 2, . . .) by u s . The spike map induces an era map u s i— > u s+ i, which recursively determines 
(u s ) se N from Uo = uq, and thereby the complete spike induced sequence (wz)/=o,i.2,... of Kasner 
epochs. 

The length of an era s is determined by the value of u s : If u s S [m,m + 1) for some m G N, 
then the length of the era is m/2, if m is even, and (m + l)/2, if m is odd. In the stochastical 
framework, in analogy with 0, let K(m) denote the probability that a randomly chosen element 
of an era sequence (u s ) s= o,i,2,... lies in the interval [m,m + 1). 

Theorem 2.1. Let (u s ) se N be a generic spike-induced sequence of eras; then the probability that a 
randomly chosen element o/(u s ) se N lies in the interval [m, m + 1) is 
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Proof. We perform the proof in three steps. 

Step 1: Representation of the era map u s H > u s +i by means of continued fractions. The spike 
map ([9]) yields the following respresentations in terms of continued fraction expansions of u± for 



the four different cases of ([9]): 

u- = [2 + fc;fei,fc 2 ,fe 3 ,...] H- it+ = [fe;fci,fc 2 ,fe 3 ,...], (11a) 

u_ = [2;fci,fc 2 ,fc 3 .fc 4 ,...] i-)- u + = [fei;fc 2 ,A;3, ...], (lib) 

U- = [1; 1,^2,^3,^4,...] !->■ u+ = [fez; fe, kl,. . .] , (11c) 

u- = [1;1 + At, fez, fc 3j — ] i-> w+ = [fe; fc 2 , £3, fei, ■ • •] , (Hd) 



where k and fcj, i = 1,2,..., are positive natural numbers. Recall that (lllb[) - (llld[) describe a 
change of era. The spike map (jTTJ) leads to a representation of the era map u s M> u s+1 : If 

u s = [x s ; k s+ i, x s+ 2j • • •] j 

then 

([x s+1 -(x s mod2);x s+ 2,x s+ 3,...] if k s+1 - (x s mod 2) ^ , 
u.s+i — \ (12) 

[[x s+2 ;x s+3 ,x s+ 4, ...] if x s+ i - (x s mod 2) = , 

when we use the standard notation (xmod2) for natural numbers (where we recall that (xmod2) = 
if x is even; (xmod 2) = 1 if x is odd). 

Step 2: Iterate ([T2"T) and represent the elements u s , s = 0,1,2,..., of the era sequence in 
terms of the partial quotients of the continued fraction representation of Uo- Let Uo be given 
by Uo = [feo; fei, fe2, • • •]• We begin by constructing from (fco, fei, fe 2 , • ■ •) a- n auxiliary sequence 
(k' , k'i, k' 2 , ■ ■ ■) according to k' = feo and 

4+1 = ks+i - « mod 2) = k s+1 - ([$^* fei] mod 2) . (13) 

In the generic case, the number will appear (infinitely often) in (k , k' 2 , ■ ■ •)■ We remove these 
entries and denote the arising (sub)sequence by (fe , k[, k' 2 , . . .). A typical example is 

(k ,h,k 2 ,...) = (4,3,6,5,1,2,1,7,2,7,1,1,3,...), 
(4, k[,k' 2 , ...) = (4, 3, 5, 4, 1, 1, 0, 7, 1, 6, 1, 0,3,...), 
(fe ,fci,fe 2 ,...) = (4,3,5,4,1,1, 7,1,6,1, 3,...). 

Accordingly, we have 

K = < +Zs = k s+Zs - ( [ Y.to' 1 h ^ m0d 2 ) ' (14) 

where z s is the number of slots (i.e., zeros) that have been removed up to index number s of the 
sequence (fe , k[, k' 2 , . ■ .). 

Comparing (fT2|) and (IT3l) we find that the era sequence reads 

uo = W a ;ki,k 2 ,---}, u 1 = [K , 1 ;k 2 ,k 3 ,...}, u 2 = [k 2 ; fc 3 , ^4, •••],•• • 

up the smallest index s for which k s+ i — (k' s mod 2) = 0; at that point, the second case of (TT2"j) 
becomes relevant, which leads to an omission of n' s+ i (= 0). The appropriate sequence that enters 
our representation of (u s ) sg N is thus the constructed (sub)sequence (fc , k[, k' 2 , . . .): Making use of 
(fc , fei, fe 2 , . . .) we are able to give the element u s of the era sequence as 

u s = W s+Za ;k s +z s +i,k s+Zs+ 2 1 . . .} = [k' g ;k s+Zs+ i,k a+Zs+2 ,...], (15) 

where z s is the number of omissions up to index number s. It follows that the probability K(m) that 
a randomly chosen element of (u s ) sg N lies in the interval [m,m + 1) is identical to the probability 
that a randomly chosen element of (fe , k[, fc 2 , . . .) is equal to m. 



3 STATISTICS 



6 



Step 3: Here we use the stochastical properties of sequences (k' Q , k[, k' 2 , . . .) that arise from 
(generic) sequences (fco, fci, hi, . . .) to derive the desired probability distribution. As a preliminary, 
we consider sums of the type X)i=o^<> w h ere (ho, k\, k%, ■ • •) are randomly distributed natural 
numbers, e.g., according to Khinchin's law ([7]). Let < p < 1 be the probability for an element 
of (fco, fci, fe, • • ■) to be odd (which is Sm=i m odd K(m) in the case of Khinchin's law). Note that 
p = 1/2 is not needed; in fact p « 0.6515 in the case of of Khinchin's law. The sum 53i=o ^» * s 
even or odd, if the number of odd summands is even or odd, respectively. Hence, making use of 
the binomial distribution, we find that the probability for 53i=o ^» ^° ^ e even or °dd is 



i= 0,2 even ^ ' i— 0,i odd 



respectively. The two probabilities are not exactly 1/2, but they converge to 1/2 rapidly as j grows, 
and hence the two probabilities are equal to 1/2 in the asymptotic limit of very long sequences H 
Finally, let Uo = [ko] k%, k%, . . .] be a generic real number. Consider the truncated sequence 
(ko, k\, . . . , k n ), where n ^> 1, and denote the probability distribution of its elements by K n . In 
the limit n — > oo, K n converges to Khinchin's distribution K, cf. ([7]). In the auxiliary sequence 
(k' , k'i, . . . , k' u ) constructed from (fco, k\, . . . , k n ), the number zero appears in approximately 21% 
of the slots. This is because of (|T3l) : k' s is if and only if k s equals 1 and J2i=i ^» i s °dd- The 
probability of the former event is K n (l) ks K(l), cf. ([7]), which is approximately 42% (because 
n » 1); the probability of the latter event is 1/2 in the limit n — > oo. Since the two events are 
independent, we obtain a probability of -£T n (l)/2 w K(l)/2 m 21% that the number appears in 
a randomly chosen slot of (k , k[, . . . , n' n ). Since the (sub)sequence (k' , k[, . . . , k' n ) is constructed 
from (kq, k' x , . . . , n' n ) by removing each appearance of the number 0, (k' , k' x , . . . , k' n ) is shorter than 
(k' , k[, . . . , K,' n ) by 21%, i.e., n w 0.79 n. We are thus able to compute the probability Ka(m) 
that a randomly chosen element k' s of the sequence (k' , k[, . . . , k' n ) is equal to m. There are two 
ways that the event k' s = m can occur. Either the corresponding element k s+Za , from which k' s is 
generated, see (fT4|) . is equal to m and X)i=o s 1 * s even ? or h 1S equal tom + 1 and X)i=o s 1 ^ 
is odd. Accordingly, 

lff n (m) + ^ n (m + l) 

** (m) = i-puii ■ (16) 

The limit n — > oo corresponds to n — > oo; since K n (m) — > K(m) in this limit, (JT)) implies 



| if (m) + i if(m + 1) ^m + 2 



#M = 1 T ,\ = log 3 —— - log, , (17) 

which concludes the proof. □ 

Corollary 2.2. Let (u s ) sg N be a generic spike induced sequence of eras; then the probability that 
a randomly chosen era in </izs sequence possesses length n is given by 

Probability (length of WE = n) =: L(n) = log 3 ( ^ + 3 ) - lo S3 + | ) ■ ( 18 ) 

Proof. This is a direct consequence of Theorem l2.ll The length of era s, where u s = [x s ; x s+ i, x s+ 2 , . . 
is x s /2 if x s is even and (x s + l)/2 if x s is odd. Therefore, for era s to be of length n, x s must be 
either 2n or 2n — 1. From (|15|) and ((TO)) we see that the probability that one of these events occur 
is K(2n) + K(2n - 1), which yields ||IS}. □ 



3 Statistics 

3.1 Comparison between BKL eras and spike eras 

Let us compare some consequences of the probability distribution ([7]), which determines the proba- 
bilities for prescribed lengths of BKL eras in BKL sequences of Kasner epochs according to eq. ©, 

3 For j = 10, the error is already of the order of 10 — 6 in the case of Khinchin's law; for a mere 100 elements, the 
error is already less than 10 -50 . 
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Sequence 


1 


2 


3 


4 


5 


10 


50 


100 


500 


era 
era 


41.50 
36.91 


16.99 
16.60 


9.31 
9.60 


5.89 
6.28 


4.06 
4.44 


1.20 
1.39 


5.5 x lCT* 
6.9 x 10" 2 


1.4 x 10" 2 
1.8 x 10" 2 


5.7 x 10~ 4 
7.2 x 10~ 4 



Table 1: Probabilities (in %) that a randomly chosen element of a BKL/spike-generated Kasner 
sequence of Kasner epochs (u s ) se rj, is in the interval [to, to + 1), m — 1, 2,3,.. .. These are the 
probability distributions K{m) and K(m), see and [fTD"]), 



Length 


1 


2 


3 


4 


5 


10 


50 


100 


500 


era 
era 


41.50 
53.50 


16.99 
15.87 


9.31 
7.75 


5.89 
4.61 


4.06 
3.06 


1.20 
0.83 


5.5 x TO" 2 
3.5 x 10" 2 


1.4 x 10" 2 
0.9 x 10" 2 


5.7 x 10" 4 
3.6 x 10" 4 



Table 2: Probabilities (in %) that a randomly chosen era/era of a BKL/spike generated sequence 
of Kasner epochs, is of a prescribed length. These are the probability distributions L(m) — K(m) 
and L(m), see © and (jl8"]l . 



and the probability distribution (|10[) . which determines the probabilities for prescribed lengths of 
eras in spike-induced sequences of Kasner epochs according to eq. ([18]), see Tables Q] and [2] As 
seen in Table HJ eras have the tendency of being shorter than BKL eras. The probability that an 
era contains merely one epoch is larger than 50%; the probability that an era consists of n > 1 
epochs is smaller than that of a BKL era. Asymptotically, for n 3> 1, we have 

Probability (length of era = n) = (log 2) -1 n~ 2 (l - 2r?r 1 + 0(n~ 2 )) , 
Probability (length of era = n) = (log3) _1 n -2 (l — n^ 1 + 0(n~ 2 )) , 

and hence the two probabilities are asymptotically proportional with a proportionality factor 
log 2/ log 3. 



3.2 Dominance of small curvature phases 

Following [1] we first introduce T > 3, where we are foremost interested in the case T>1. Then 
we define a small curvature phase of a BKL or spike induced sequence of Kasner epochs (mj)z g n as a 
connected and inextendible piece £cfJ such that it; > T V/ G C. During a small curvature phase 
ui is thus monotonically decreasing from a maximal value by the BKL map to a minimal value in 
the interval (T,T + 1], while the spike map yields a monotonic decrease from a maximum value 
to a minimal value in the interval (T, T + 2]. The complement of the concept of a small curvature 
phase is a large curvature phase, which is defined as an (inextendible) piece of the sequence of 
Kasner epochs such that it/ < T for all I. 

During small/large curvature phases, the 'Hubble-normalized' curvature is comparatively small/large, 
see eq. (|25p below. While a small curvature phase can be viewed as an era/era that is terminated 
prematurely at T, a large curvature phase typically consists of many eras/eras; clearly, small and 
large curvature phases occur alternately. In the following BKL example, where the choice T = 4 
has been made, the large curvature phase contains two and a half eras. 

small curvature phase large curvature phase 

1.14 -)■ 7.29 -»■ 6.29 -»■ 5.29 -»■ 4.29 -»■ 3.29 -»■ 2.29 -»■ 1.29 -»■ 3.45 -»■ 2.45 -»■ 1.45 -»■ 2.24 -»■ 1.24 



Combining the probabilistic viewpoint with the concept of small/large curvature phases lead 
to a fundamental result in the description of the BKL and spike induced Kasner sequences. For 
generic Kasner sequences (u/)zgN> small curvature phases dominate over large curvature phases in 
the following sense: 

Theorem 3.1. Let oe o, generic BKL or spike-induced Kasner sequence and let T be 

arbitrarily large. Then for a randomly chosen epoch u the probability for the event u > T is one 
and the probability for the event u < T is zero. 
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Proof. Let us first consider the BKL case and let us for simplicity choose T G N. We let uq = 
[fco; fci, fe) • • •] be a well approximable irrational number and consider the BKL Kasner sequence 
of epochs (w;); £ n and the associated era sequence (u s ) s£ n (where we recall u s — [k s ; k s+1 , ...]). 
Consider the truncated sequence of eras (u s ) s <„ and denote the associated probability distribution 
by K n (which converges to Khinchin's law K in the limit n — > oo). Let u be a randomly chosen 

element from the associated (finite) sequence of epochs (u/);<fe H \-k n - We denote by £ m the 

collection of epoch of eras of length to and by P n (u&£ m ) the probability that the epoch u belongs 
to an era of length to; it is given by 

P n (ue£ m ) = (J2 m ' Knim^y 1 mK n (m) . (19) 

m' 

The probability for u < Y, which we denote by P n {u<T), is thus obtained from 

T— 1 oo T — 1 

P n (u<T) = J2 Pn(ue£ m ) + J2 ^»Mm) — 



m=l m=T 



Y,m'K n {m') 



T-l oo 

J2 mK n {m) + (T - 1) £ K n (m) 

m=l rn=~T 



For a generic number uq, in the asymptotic limit n — > oo, the probability K n (m) converges to 
K(m) given by Khinchin's law ([7]). Since 

K(m) = — - [ m - 2 - 2m- 3 + 0(m" 4 )] (m -» oo) , (20) 

the sum X)m=r ^n( TO ) approaches the convergent series X)m=r K( m ) m the limit n — > oo; however, 
the sum ^ m mK n (m) diverges as n —> oo, so that 

P(u<T) = lim P n (u<T) = 0, (21) 

v ' n— »oc v ' 

irrespective of the value of T. We conclude that the probability that an epoch u, which is randomly 
chosen from the Kasner sequence (itz)j e N, belongs to a large curvature phase (u < T) is zero; with 
probability 1 a randomly chosen element belongs to a small curvature phase (u > T). 

It remains to consider the recurring spike case, where we choose an arbitrarily large T such 
that T = 2v + 1 with v G N for simplicity. Let L n (m) denote the probability for an era of the 
truncated sequence to be of length m; L n (m) converges to L(m) as n —J> oo, cf. (fT8]) . Then 

P n {ue£ m ) = (^2,m! L n {m')\ minim), (22) 

m' 

and the probability for u < T reads 

p ( V* p t c e \ . V" p r c p ^ u Sm^i TO -^n (to) + u E^+i L n (m) 

P n lu<T) = ) P n (u€£ m ) + } P n {ue£ m )— = ^ , f , A ■ 

* — ' * — ' to > . m L n (m ) 

m=l m=v+l L -' m ny 1 

In the case of a generic sequence, L n (m) converges to L(m) as n — > oo, where 

L(m) = - - [to -2 - to ^ + 0(toT 4 )] (to -> oo) . (23) 

log o 

Hence the sum X)m=u+i -^n( m ) converges as n — » oo. However, the sum J2 m mL n (m) diverges as 
n — > oo and P(u<T) = lhxin^oo P n (u<T) =0. □ 

The underlying reason for the dominance of small curvature phases is hence the failure of the 
probability distributions 0, (fTO]) (and jSJ, (fT8|l). to generate finite expectation values, since 



mK{m) — m L(m) = oo , mK(m) = oo = 

m— 1 m—1 m—1 m—1 



which is due to the infinite tail of the distributions. Accordingly, the average length of an era/era 
is ill-defined. 
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3.3 Weyl stochastics 

The magnetic part of the Weyl tensor for any Kasner state is zero, and hence the Hubble- normalized 
Weyl scalar is given by = ^H~ 4 E a pE a P , where E a p is the electric Weyl tensor and H is the 
Hubble variable, which is related to the expansion 9 of the timelike reference congruence according 
to H = (9/3, see [4J [TO] . Expressed in terms of the Kasner parameter u, is given by 



W£ = -81 PlP2P3 = 1 2 7 . 25 



The behavior of over generic sequences («z); e N of Kasner epochs is statistically described by 
the following Theorem: 

Theorem 3.2. Let U\ = U\ = [k±; fe,^, . . .] be a generic real number, (u;)igN the associated BKL 
or spike induced sequence of epochs, and (u s ) se N the associated BKL or spike induced sequence of 
eras. Consider the truncated sequence (ui,...,u„) of eras/eras. Then the average (Wj|)„ of Wj| 
over all epochs of the first n BKL and spike-induced eras satisfies 

(W 2 )„ = 19.580317157(8) - J— , (26) 

log 2 n l + x(n) 

where \ depends on if we consider an era or and era, but in both cases it is a function satisfying 
liminfn^oo x( n ) — an d nmsu Pn^oo xi n ) = 00 • 

Remark. Although the results are similarly the same for BKL and spike-induced sequences of 
Kasner eras, note that the concepts of eras differ for the two cases. 

Remark. The function \ does not converge to zero exactly, but 'almost'. Specifically, x( n ) converges 
to zero as Af 3 n — >• oo, where Af is a subset of N minus a subset of logarithmic density zero. The 
average (W£)„ (slowly) converges to zero as n —> oo. Intuitively, this is due to the dominance of 
small curvature phases over large curvature phases. 



Proof. We first prove Theorem 13.21 for the BKL case and then for the recurring spike case. In the 
BKL case we have: 

Lemma 3.3. Let U\ = u\ = [ki] k 2 , k%, . . .] be a generic real number, (m;); 6 n the associated BKL 
sequence of epochs, and (u s ) s£ n the associated BKL sequence of eras. Consider the truncated 
sequence (ui, . . . , u„) of eras. Then the average (Wj|}n of over all epochs of the first n eras 
satisfies 

(Wl >„ - 19.580317157(8) — L- X + ° (1) , (27) 

log 2 n 1 + xN 

where % zs a function satisfying liminfn^oo x( n ) — an d umsu Pn->oo x( n ) = °°- 

Let n ^> 1. The average of over all epochs (= k\ + ■ ■ ■ + k n ) of the first n eras is 

n feiH h^n n n k s 

(^) n = (j2^y 1 E w 2 w = (E fc _1 EE w K(K}+ J ), 

i=l 1=1 i=l s=l j—1 

where {u s } = [0; fc s +i, . . . ] = is the fractional part of u s ; recall that the integer part [u s ] is 
equal to k s . Since n > 1 we may consider {u s } =: a as a random variable on (0, 1); its probability 
distribution is w(x) = (log2)~ 1 (l + x)' 1 , see [3] with x = u _1 and w(x)dx = —p(u)du. We are 
thus able to replace YljLi Wk({ u s} +i) by its weighted average W^(k s ) which is given by 



W 2 (m) := / w{x) Y\ W&(j + x)d>c = / w(x)W£(j + x)dx . 
Jo j 

Accordingly we obtain 



m „i 

]Wi{j + x)dx = Y, / ' 

1 .7=1 J ° 



<W|>n = (t^rtWiks) = ^"pifW , (28) 

i=l »=i n 2^1=1 K i 
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where K n denotes the probability distribution associated with (k\, . . . , fc„), which converges to K 
as n — > oo. 

To compute the numerator of (|28p we first note that 



l 

2 , 



w(x)Wk(j + = / io(x)W|(j + x)d>r + / w(x) [W|(j + x) - Wk(j + x)] dx 



+ x) - 162j- 3 / io(x)(x - x)dx + 0(j" 4 ) 

J + Sr+1/2 

W2(x)dx + 0(r 4 ) (29) 

j + Sr-1/2 

as j — >■ oo, when we use X = J Q xw(>c)du = (log2) _1 — 1. Accordingly, for m > 1/ > 1 we find 
W£(m) - W|(M) = ^ / w(x)W£(j + x)cte = / + 0{M- 3 - m~ 3 ) 

j=M+l^ Q JM+a+1/2 

81 

= SlfM- 1 - m _1 ) - — — (M- 2 - m- 2 ) + 0(M- 3 - m- 3 ) . 
v ' log 2 v ' 

Since J2 m K n (m)W^(m) converges to J2 m K(m)yV^ i (m) as n — » oo, we consider the latter series. 
Choosing N 3 M » 1 we obtain 

00 M 00 . Q1 01 1 

^if(m)W^(m) = ^A:(m)W|(m)+ ]T tf(m) (vt£ (M) + - - - + 

m=l m=l m =M+l 

m— 1 

A numerical calculation then yields the result 19.580317157(8). 

It remains to analyze the denominator of (128[) . which is the mean length of the first n eras, or, 
in other words, the arithmetic mean of the first n partial quotients of [k\; k-z, ■ ■ ■ ]. From the theory 
of continued fractions we have the result that 

n 

J2ki=n log 2 n (1 + X (n)) , (30) 
»=i 

where x is as described in the lemma and the remark; see, e.g., [H1[T3]. By combining (|30p with 
the result of the previous calculation the lemma is proved. 

Let us now turn to the recurring spike case. By considering the union of all Kasner sequences 
and noting that the spike Kasner map is just the square of the BKL map it is perhaps intuitive 
that we will obtain a similar result as in the BKL case, but let us give a formal proof nonetheless. 

Lemma 3.4. Let u± = u± = [ki;k2,ks, . . .] be a generic real number, (it;);gN the associated spike 
induced sequence of epochs, and (u s ) s£ n the associated spike induced sequence of eras. Consider 
the truncated sequence (ui, . . . , u„) of Was. Then the average (W^) n o/ W£ over all epochs of the 
first n eras satisfies 

(W&>„ - 19.580317157(8) — , (31) 

log 2 n 1 + XW 

where x is a function with the same properties as \ of Lemma \3.3\ . 

Remark. The statement of this lemma is formally identical to that of Lemma [3~31 which motivates 
the formulation of Theorem 13.21 Note, however, that the concept of an era of a spike induced 
sequence of Kasner epochs differs from that of a BKL sequence. 

For the proof of the above lemma we need P^)) - p31) ; in particular we recall that u s = 
[k' s , fc s+Zs+ i, . . . ]. Let us denote the length of era s by l s , where l s = l(k' s ) with 



l(k') = { % k ' 6Ven ' (32) 
K J 1 iL+I y odd , 
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The average of W£ over all epochs (= l\ + • • • + l n ) of the first n eras is 

n \-l n n n l s 

(Wl)„ = (^Z i )- 1 Yl W^u l ) = (Y / kr 1 Y,Y, W K&-( k >od2) + {u s }). 

i=l 2 = 1 i=l s=l j = l 

In complete analogy with the proof of Lemma 13.31 we are able to replace the inner sum by by the 
weighted average W^fcg), where 

W K( fc : = / Y] w k( 2 J - (k' mod 2) + >c)d>c = V / w(x)W^ (2j - (k' mod 2) + x)dx . 

Accordingly, 

0«). -(£<,)- t^W)- ^^?^ . W 
»=1 s=l n <oi=l H 

where -K„ denotes the probability distribution associated with (fci, . . . , k n ), which converges to K 
as n — > oo. 

To compute the numerator of ([33| we use (129|) . where j is replaced by 2j — fc' mod 2. In addition 
we make use of the relation 

2j+a+l/2 -, .2J+0-1/2 i /-2j+o+l/2 i ^2j+a+3/2 

Wl{x)dx = - w£dx + - widx + - W&x + 0(r 4 ), 



2j+a-l/2 4 J2j+a-3/2 z J2j+a-l/2 4 J2j+a+l/2 

which holds for all a (e.g., x - fc' mod 2). Let A;' > K' » 1 with fc' mod 2 = K' mod 2. We find 

i(fe') 



W£(fc') - WKif') = J2 / w (^) W k(2j- (fc'mod2) + x)dx = 



j=J(^')+l 

= t / WK^dx + - / Wl{x)dx + - / W> (x)dx + o{(i<t 3 ) 

= - (kT 1 ) - ( 2 + lo S 2 ) ((^')- 2 - (fc'r 2 ) + 0((^')- 3 ) 

Since J2k> ^n(k')yV^(k r ) converges to J2 k , K(k')W^(k') as n — )• oo, we consider the latter series. 
Choosing N 9 ff ' > 1 we obtain 



£ ^(fc')Wl(fcO = £ ^( fe o(wl(^-i) + |(-^-i)+o(^- 

:' odd 



fe'=K'+l k'=K'+l 

k' odd 



fc' even 

and thus 

2 K(k')Wl(k') = £ - Wl(K< - l) l 0g3 jL±± - Wl{K>) log 3 jL±| 

fc'=l k' = l 



— —\0 K ' + 1 81 1 Qf_}_ 

2 K> 0&3 K> + 3 21og3 K 12 + ^K> 



A numerical calculation then yields the result 12.35380467921(1). 

It remains to analyze the denominator of (|28[) . which is the mean length of the first n eras. 
Using (|3"2"1) we have 
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where log 3 2 is the probability for a randomly chosen element of (fcj, k' 2 , ■ . . ) to be odd. From the 
construction of the auxiliary sequence k' 2 , . . . , «4) and the sequence (fc 1; k' 2 , . . . , k' n ), see (fT3|) - 
(|15p . we recall that (fc^, fe^, - - - , k' n ) is shorter than k' 2 , ■ ■ ■ , ) by the fraction ± A'(l) = 1 - 
(21og 3 2) _1 ; accordingly, we have n = 2(log 3 2) n. Furthermore, we obtain 

i—1 i—1 i—1 



because, on the average, the elements of (/Cj, . . . , k' a ) are smaller than the elements of (fci, . . . ,k„ 
by one half; we again refer to the proof of Lemma |2~T1 Therefore, 



n 1 n 11 1™ 1 

5Z Zl= 2^-' fc4 '4" + 2 n l0g3 2= 2^ fcl= 2" l0g2 " ^ + X< ™^ (34) 
= n log 3 2(1- log 2 log 2 3 + log 2 n) (l + x(n)) = n log 3 n (l + x(ra)) 

Noting that 12.35380467921(1) log3 (log2)- 1 = 19.580317157(8), the lemma is proved. This con- 
cludes the proof of Theorem 13.21 □ 



Remark. The results of Lemma 13.31 and 13.41 differ slightly when expressed in terms of a different 
parameter; instead of n we use the number N of epochs up to (and including) era/era n. Let us 
restrict ourselves to the main point. In the context of Lemma l3.3l we have N = X)"=i ~ n 1°§2 n i 
cf. (|3U| . Therefore, log 2 n = (log 2)~ 1 W(N log 2), where W is the Lambert W function (product 
logarithm). For N ^> 1 we thus have 

log 2 n = (log 2)- 1 ( log(iV log 2) - log \og(N log 2) + o(l)) 
= (log 2)- 1 ( log N -log log TV + log log 2 + o(l)) . 

Analogously, in the context of Lemma [3.41 we have N = h — n log 3 n, cf. (|3"4"|) . Therefore, 

log 2 n = (log2)- 1 lF(A^log3) and 

log 2 n = (log2) _1 (logA r -loglog7V + loglog3 + o(l)) . 



4 Concluding remarks 

The BKL and the BKL era maps have attracted a lot of attention in the literature, perhaps 
due to the tantalizing chaotic features that they indicate, see [TT] and [H]-[3n], an d references 
therein. The present paper opens up similar exploration possibilities for recurring spike formation 
— we have only given a few examples of what can be done. As in the BKL case, one can derive 
probabilities for other quantities, but one can also use the explicit solutions that describe spike 
transitions to attempt to extend the associated maps in perhaps similar ways as have been done 
for the BKL case. Another possibility is to focus on the state space picture in [9l[T0], where [9] give 
some examples of so-called concatenated chains, where axes permutations are not quoted out (for 
a discussion, see |10j). In particular, one can explore periodic spike chains, which are the recurring 
spike analogues of the so-called heteroclinic cycles in the BKL case, which recently have attracted 
attention [3TII32]. 
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